Generalized PSK in Space Time Coding 



Guangyue Han 

Department of Mathematics 
University of Notre Dame 
Notre Dame, IN 46556. 

e-mail: Han.13@nd.edu 
URL: http:/ /www. nd.edu/~eecoding/ 

January 14, 2003 



Abstract 

A wireless communication system using multiple antennas promises reliable trans- 
mission under Rayleigh flat fading assumptions. Design criteria and practical schemes 
have been presented for both coherent and non-coherent communication channels. In 
this paper we generalize one dimensional phase shift keying (PSK) signals and introduce 
space time constellations from generalized phase shift keying (GPSK) signals based on 
the complex and real orthogonal designs. The resulting space time constellations real- 
locate the energy for each transmitting antenna and feature good diversity products, 
consequently their performances are better than some of the existing comparable codes. 
Moreover since the maximum likelihood (ML) decoding of our proposed codes can be 
decomposed to one dimensional PSK signal demodulation, the ML decoding of our 
codes can be implemented in a very efficient way. 

Index Terms-space-time coding, multiple antennas, orthogonal designs, diversity, phase 
shift keying 



1 Introduction and Model 



Multiple antenna communication systems have been actively studied recently By 
exploiting the temporal and spatial diversity at both transmitter side and receiver side, such 
kind of systems can increase the channel capacity compared to single antenna communication 
systems, consequently promising more reliable data transmission for high rate applications. 

We investigate a communication system under Rayleigh flat fading assumptions and we 
further assume fading is quasi-static over a time period of length T. Denote by M or N 
the number of transmitting or receiving antennas, respectively. Let p represent the expected 
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signal-to-noise ratio (SNR) at each receiving antenna. For the system above, the basic 
equation between the received signal R, which is a T x N matrix, and the transmitted 
signal $, which is chosen from a T x M matrix constellation V = {$1, $2, • • • , $l} (L is the 
constellation size), is given through: 



_P_ 

M 



R = JJLqh + W, 



where the M x N matrix H accounts for the multiplicative complex Gaussian fading co- 
efficients and the T x N matrix W accounts for the additive white Gaussian noise. The 
entries h m ^ n of the matrix H as well as the entries w t>n of the matrix W are assumed to have 
a statistically independent complex normal distribution CAf(0, 1). One can verify that the 
transmission rate is determined by L and T: 

R _ log 2 (£) 
T 

When the fading coefficients are unknown to the transmitter however known to the 
receiver, it is proven [21 EI] that the channel capacity will increase linearly with min{M, iV}. 
Since H is known to the receiver, for a received signal R the maximum likelihood (ML) 
decoder will take the following evaluation to resolve the most likely sent signal (codeword): 



$ = arg min 



M 



where || • || represents Frobenius norm. Let P$ u $ , denote the probability that ML decoder 
is mistaking $/ for The upper bound of this error probability has been derived in |Hj: 



j M r • , ■ - - Y 

p ^i> - 2 n 

771=1 



AM 



Through the analysis of the above upper bound, design criteria for designing codes for 
the coherent channel with ideal channel state information (CSI) are proposed in |Uj: 

The Rank Criterion: In order to achieve the maximum diversity MN, the matrix B(<&i, $J) 
$z — $;/ has to have full rank for any codewords and <ly. If B(<&i, $//) has minimum rank 
r over the set of two tuples of distinct codewords, then a diversity rM is achieved. 

The Determinant Criterion: Suppose that a diversity benefit of rM is our target. The 
minimum of rth roots of the sum of determinants of all r xr principal cofactors of A(&i, $;/) = 
B(<&i, <&ii)B*(<&i, $//) taken over all pairs of distinct codewords and <Ev corresponds to the 
coding advantage, where r is the rank of A($i, $//). 

Note that for the case when T = M and full diversity is achieved, the above design 
criteria can be simplified as follows: 

Construct a constellation of matrix V = {$1, $2, • • • , $l} such that the diversity prod- 
uct 

nV = mmi|det($,-$,Or /M 

is as large as possible. 
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Orthogonal designs have been investigated for the constellation construction for coherent 
channels. Using a complex orthogonal design, Alamouti [Tj proposed a very simple trans- 
mitter diversity scheme with 2 transmitting antennas. A very interesting coding scheme 
from the real orthogonal designs is presented in [Hj. Applying the similar idea as in pQ, the 
authors also explain how the symmetric structure of the orthogonal design codes leads to a 
much simpler decoding algorithm. However no explicit constellations have been studied in 
details in these work. 

For a slow fading channel, for instance a fixed wireless communication system, the coher- 
ent assumption is reasonable because the transmitter can send pilot signals which allow the 
receiver to estimate the fading coefficients accurately. However in certain situations, due to 
the difficulty of measuring the fading coefficients under the practical communication envi- 
ronment (for instance the limited resources or a fast fading link) the above coherent model 
would be questionable. Because of this, researchers put their efforts into the case when the 
channel is non-coherent, i.e., the CSI is not known by either the transmitter or the receiver. 

In jHj, Hochwald and Marzetta study unitary space-time modulation for a non-coherent 
channel. We will use the same notations as in the coherent case, so the basic equation will 
be still the same: 

R = . JL&H + W, 
V M 

however in non-coherent scenarios it is assumed that the receiver does not know the exact 
values of the entries of H (other than their statistical distribution). Another difference is 
that the signal constellation V := {$1, $2, • • • , $l} has unitary constraints: &* k &k = Im f° r 
k — 1, 2, • • • , L. The last equation simply states that the columns of form a "unitary 
frame", i.e., the column vectors all have unit length in the complex vector space C T and the 
vectors are pairwise orthogonal. The scaled matrices y/T&k, represent the codewords used 
during the transmission. 

The decoding task asks for the computation of the most likely sent codeword $ given 
the received signal R. Under the assumption of above model the ML decoder will have to 
compute: 

$ = arg max H-R*^!! 
<i>e{*i,<i>2,-- 

for each received signal R (see jH]). 

It has been shown in jS] that the pairwise probability of mistaking for $// using ML 
decoding satisfies: 



1 n 

< 2 n 

m=l 



1 + MMiiW"' v 

4(1 + pT/M) 



where 5 m (^^i>) is the m-th singular value of An important special case occurs when 

T = 2M. In this situation it is customary to represent the unitary matrix $^ in the form: 

Note that by the definition of <3>fc the matrix $^ is a M x M unitary matrix. In most of 
the literature mentioned above researchers focus their attention on constellations having the 
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special form A differential modulation scheme is discussed in [3] and this special form 

is used. In this scheme, one does not send the identity matrix (the upper part of the signal) 
every time. Instead of sending 

f hi\ ( hi \ ( hi \ 

one sends 

This increases the transmission rate by a factor of 2 to: 

= log 2 (L) = 2 log 2 (^) 
M T 

Let \I/ T and R T denote the sent and received signals at time r respectively, the ML decoder 
of the above differential space time modulation scheme will have to compute: 

^f T = argmin ||i? T — \E'i? r _ 1 ||. 
For the constellations with the special form Ijl.ljl . the pairwise error probability satisfies: 

1 n 

m=l 

As explained in [Sj, at high SNR scenarios the right side of the above inequality is governed 
by the diversity product: 

nV = rnini|det(vI/ i -vI/ r )| 1/M . 

So the unitary differential modulation design criterion for a non-coherent channel is to choose 
a constellation V such that J} V is as large as possible. We call a constellation V fully diverse 

if nv > 0. 

Interestingly enough the normalized complex orthogonal design Alamouti codes and real 
orthogonal design codes can be used in non-coherent scenarios too. Our work will be mainly 
about non-coherent channels, however without any doubts the resulting codes can be used 
for coherent channels too. In this paper we will show how to construct space time codes 
from these schemes using generalized phase shift keying (GPSK) signals. Our results can be 
applied to generalized orthogonal design in jjj. 

This paper is organized as follows: in Section |21 three series of 2 dimensional GPSK 
constellations from complex orthogonal designs will be presented. An algebraic calculation 
will show that they have larger diversity products than the original orthogonal design con- 
stellations. As a consequence of larger diversity, the ML decoding of a GPSK constellation 
gives better performance. In Section 3, we explicitly construct GPSK constellations from 
the real orthogonal designs. Fast decoding algorithms for the proposed constellations will be 
presented in Section 2 and Section 3. Finally in Section 4 we will give the conclusions and 
some future work. 



1 + 



4(1 + 2p) 



1.21 
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2 GPSK Constellations from the Complex Orthogonal 
Design 

A very simple yet interesting complex orthogonal scheme is described in pQ. We will consider 
the normalized version of this proposed code, i.e., every element of this code is a matrix given 
by 

where \a\ 2 = \b\ 2 = 1. Observe that 0(a, h) is a unitary matrix and a constellation 0(n) 
with size L = n 2 is obtained by letting a and b range over the set of n-th roots of unity 



= {l,e 



2m /n 



2ni(n—l) /n 



namely the entries of the matrix are chosen from scaled one dimensional PSK signal set t?. 
This would be the most commonly implemented Alamouti's scheme. The diversity product 
of the constellation 0(n) is 

Note that 0{n) is similar to d in the sense that all the elements in 0(n) have unit energy 
and every pair of elements differ only by the phases. 0(n) is a subset of the special unitary 
group 



SU(2) 









{(-!• 




\a\ 



\bf 



In this section we will present three series of unitary constellations as finite subsets of SU(2). 

The basic principle to decode the constellation 0(n) has been discussed in [10J. We are 
going to describe this decoding process in another way with more details and generalize it 
to our constellations. Consider a non-coherent wireless communication system modulated 
by 0(n) with 2 transmitting antennas and N receiving antennas and assume the differential 
unitary space time modulation |H] is used. Let X and Y denote the received matrices at time 
block r and r + 1, respectively, then the ML decoder will perform the following decoding 
task 



(a, 6) 



arg mm 



a 

-b* 



X 



With simple matrix manipulations, one can check that 



N 

(a, b) = arg min > 

»=i 



Y 2i 



J_ / a b 
V2 V ~ b * a * 



Xu 
X21 



N 

arg min > 

1=1 



arg mm 



N 



a,b&§^-^ \X 



It 



2; 



2/ 



1 
71 



1 / Xij A 2 j 



Xu 
v* 

A 2i 



Xu 



2i 



1 2 



a 
6 



5 



■ f 2 

arg mm > * ^ 

-^tt \X u \ 2 + \X 2i \ 2 



1 1 1 2 2 



+ 



1 1 1 2 2 

Further algebraic simplifications show that ML decoding is very simple: 



iV 



N 



, b) = arg max 



(2.1) 



i=i 



i=l 



where ^ = A^Y^ + A^Y^* an d Wi — X^Yu + AijY 2 *. Since a and b are independent of each 
other (a and 6 can be chosen freely in the set the evaluations above amount to 

rt 

a = arg maxi?e(a Zj), 



and 



b = arg max Re(b^~^^ Wi). 



i=i 



Rewrite a as a = e 27rjfl//n and b as b = e 2nk ^ n and let 

Lrl = Lr + 1/2J, 

i.e., [ r l denotes the smaller of the closest integers to r. The ML decoder will take the 
following simple form: 



(2.2) 



nargZl - nargW 

, k = 

2tt ' [ 2tt 

where Z = J2ti Z i and W = Eti W i- 

Assume that a communication channel is modulated by one dimensional PSK signals 
Let a = e 27r '/ ra £ $ denote the sent signal at one time and Z denote the corresponding 

corrupted signal. The ML decoding will look for the closest signal d = e 27T ^ n to a in the 

signal set i.e., 

n axgZ 



J 



2tt 



(2.3) 



Compare Formula ()2.2|) and Formula (|2.3|) . one can conclude that the decoding of space 
time code 0(n) is decomposable and can be implemented by taking one dimensional PSK 
demodulation twice. The normalized Alamouti codes 0(n)'s admit very simple decoding 
algorithms, since roughly only 4N complex multiplications and 4N complex additions are 
needed. We are going to introduce other complex orthogonal codes using generalized PSK 
(GPSK) signals in the sequel. Instead of allocating the same energy to every transmitting 
antenna, we attempt to exploit the transmit diversity and optimize the power allocation. 
Compared to 0(n), these codes have larger diversity, which promise good performances with 
ML decoding. Due to the symmetry from GPSK signals, these codes are also decomposable, 
therefore fast decoding algorithms can be also applied to these codes. In the following three 
subsections we present the constructions of these three series of constellations. 
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2.1 Construction 1 

Let n be an even number and let < r < ^ be the root of the following equation 

^-r) +(y^- VT^] = 4r 2 sin 2 — . (2.4) 

Consider the following sets of the scaled one dimensional PSK signals: 

Ai(ra) = {v / 2/2e i ^ fc |fc = 0, 1, • • • ,n-l}, 

A 2 (n) = {r e ^ fc |fc = 0,l,...,^-l}, 

Aa(n) = {VT^e l ^ k \k = 0, 1, • • • , n - l} . 
Consider the following subsets of 577(2): 

CiW = |(^ ^ |ae Ai(n),6e Ai(n)|, 

C2(n) = |^_^ ^ |aG A 2 (n),foe 4 3 (ra)j, 
C7 s (n) = {(_? 5 JJlae^iG^n)}. 

Theorem 2.1. Let 

Vi(n):=Ci(n)|JC2(n)UC3(n), 
Vi(n) is a fully diverse constellation of 2n 2 elements with diversity product: 



Yl Vi(n) = min <j ^ sin ^, I*/ (2^ _ r )2 + (^ - Vl - r 2 ) 2 



> . 



Proof. Take two elements A 6 Cj and B <E Cj. Without loss of generality, we further assume 
i < j. If AS G Ci(n), then 



// n 



|det(A-S)| > (v^sin-) 2 = 2 sin 
If A G Ci(ra), S G C 2 (n) or A G Ci(n), 5 G C 3 (n) then 

| det(A - 5)| > (— - rf + (— - v / r^ 2 ) 2 = 2 - v^r - V2~VT^. 
If A G C 2 (n) and B G C 3 (n) then 



\det{A-B)\ > 2(Vl -r 2 



r . 
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If A, B G C 3 (n) or A, B G C 2 (n) then 

2-7T 

|det(A-S)| > (2rsin— f. 

n 

Use the fact that r is the root of Equation (J2.4)) and compare the lower bounds of the diversity 
product in all the above cases, the claim in the theorem can be established. □ 

The following table shows how the diversity product of V\{n) compares with the diversity 
product of 0(n) when n < 12. 



n 


r 


f sinS = nO(n) 




nVi(n) 


4 


0.259 


0.5 


0.259 


0.259 


6 


0.284 


0.353 


0.246 


0.246 


8 


0.321 


0.271 


0.227 


0.227 


10 


0.360 


0.219 


0.209 


0.209 


12 


0.386 


0.183 


0.193 


0.183 



Corollary 2.2. For n > 12, 



Proof. For n > 12, 



T7 Vl(n ) = ^r sin -- 

u 2 n 



1 //v2 NO /v2 > 27T 7T 7T 7T 7T 

-A / ( r) 2 + ( vl — t)=t sm — = 2r sin — cos — > 2 x 0.386 cos — sin — 

2V 2 ;v 2 ; n n n~ 12 n 



^V2.7T 

> sm — . 

~ 2 n 



Consequently we have 



m, t \ .\V2.nl V2 V2 f I y/2 . vr 

HVi(n) =mm< — sm -, -\l (— - r) 2 + (— - VI - r 2 ) 2 V = ysm-. 

□ 

The above corollary indicates that for n > 12, the GPSK constellation V\{n) has the 
same diversity product as the orthogonal constellation 0(n), while it has twice as many 
elements: V\{n) has 2n 2 elements whereas 0(n) has n 2 elements. 

Similar to the case for 0(n), the ML decoding of Vi(n) boils down to Formula ()2.1|) . 
However we can not separate the estimation of a, b using the simple Formula (|2.2|) . because 
generally a and b are not independent anymore. If we restrict the decoding evaluation in a 
particular Cj(n), then arg(a) and arg(6) can be chosen freely. Namely within the restricted 
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searching area Ci(n), a and b are independent of each other. Assume the originally sent 
codeword falls in Ci(n), one can use Formula (|2.2j) to resolve a candidate (ai,&i): 

n arg W 

Similarly we can have candidates (ctj, 6j) for Cj(n), where z = 2, 3. The final ML decoder will 
resolve the most likely sent codeword: 

(a,b) = arg max (Re(aZ) + Re(bW)) . 

Again we conclude that the decoding of Vi(n) is decomposable and therefore V\{n) admits 
a simple decoding. The above evaluations require roughly AN complex multiplications and 
AN complex additions, therefore the decoding of the GPSK constellation Vi(n) has the same 
complexity as that of 0(n). 



nsngZ 
2tt 



2.2 Construction 2 

Let n = 2m and consider the following sets consisting of scaled one dimensional PSK signals: 

A 1 (n) = (re i % i \j = 0,l r -- ,m-l\, 



A 2 (n) = <Vl- r 2 e i{ ™ j+ ™ ] \] = 0, 1, ■ ■ • , m - l) , 



A 3 (n) = |Vl -r 2 e L ™ j \j = 0,1,--- ,m- 1 
A 4 (n) = |re i( |j = 0, 1, ■ • • , m - l| 



where 



'2sin 2 ^ + 2 v / 2sin-^ + 2 
Based on the above signal sets, construct the following subsets of SU(2): 

C ^ n ) = {(-& \ ^\ a ^ Mn),be A 2 (n)^, 

C 2 (n) = { C " g * ) |o e A 2 (n), 6 G A 1 (n] 

i b 
-b a 

i b 
-b a 



G M = {[ \ t )\ae A 3 (n),be A 4 (n] 



C ^ n ) = {[ \ t )\ae A A {n),be A 3 (n] 
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Theorem 2.3. Let 



V 2 (n) ■^C 1 {n){jC 2 {n){]C z {n)\JC i {n), 
V2(n) is a fully diverse constellation of n 2 elements with diversity product 



TT "^(n) = mm 1 r sm — ' sm — 1 ■ 

- LJ - [ n n j 



Proof Take two elements A G Cj(ra) and B G Cj(n). Without loss of generality, we can 
further assume i < j. If 2 = j, we have 



ldetM-5)] > |re^- re *f( fc + 1 )| 2 = 4r 2 sin 2 -. 

m 

If A G Ci(n), B G C 2 {n) or A G C 3 (n), B G C 4 (n), one can check that 

I det(A -B)\> 2\re l ^ k - VT^e^+^l = 2(1 - 2v / T 3 ^ 2 r cos — ). 

m 

Similarly if A G d(n), B G C 3 (n) or A G C 2 (n), 5 G C 4 (n), 

|det(A-S)| > 2(Vl - r 2 - r) 2 = 2(1 - 2y/l - r 2 r). 
If A G Ci(n), 5 G C 4 (n) or A G C 2 (n), 5 G C 3 (n), 

I det(A - 5)| > |re^ fc - re^ + ™| 2 + \Vl - r 2 e l ™ k - y/l - r 2 e i '^ + ™\ 2 = 4sin 2 — . 

2m 



It follows from the definition of r that 

2(1 - 2Vl - r 2 r) = 4r 

and naturally we will have 

7T 



2 • 2 

sin 



7T 



m 



2(1 - 2Vl - r 2 r cos — ) > 2(1 - 2^1 -r 2 r). 
m 

Compare the lower bounds in all the cases and take the minimum of them, we establish the 
claim in the theorem. □ 

The following table shows how the diversity product of V2(n) compares with that of 0(n) 
when n < 14. 



n 


r 


r sin — 


7T 

sm - 

n 


nv 2 (n) 


UO(n) 


4 


0.383 


0.383 


0.707 


0.383 


0.5 


6 


0.410 


0.355 


0.500 


0.355 


0.354 


8 


0.447 


0.316 


0.355 


0.316 


0.271 


10 


0.479 


0.282 


0.309 


0.282 


0.219 


12 


0.505 


0.253 


0.259 


0.253 


0.183 


14 


0.527 


0.229 


0.222 


0.222 


0.157 
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Corollary 2.4. For n > 14, 

IT V 2 (n) = sin-. 
±x n 

Proof. According to the definition of r, 

2(1 - 2Vl -r 2 r) = 4r 2 sin 2 — . 

m 

One can also check that 



4r sin — = 16r sin — cos = 4r cos (4 sin — ) 

m 2m 2m 2m 2m 

4 COS 2 . „ 7T , 4 COS 2 ,,.,)[. „ . o 7T 

-(4 sin 2 — ) > — — (4 sin 2 — ) > 4 sin 2 



2 sin 2 z- + 2y/2sm- + 2 ' 2m y 2 sin 2 f + 2^ sin £ + 2 ' 2m y ' 2m' 

Consequently one concludes that 

IT V2 (n) = min { r sin — , sin \ = sin — . 

- LJ - I m 2m J n 

□ 

For n > 14, the constellation V2(n) has as many elements as the constellation 0(ri), 
however its diversity prodcut is larger than that of 0(n) by a factor of a/2, i.e., n^2(^) = 
\/2 f| O(n). The corollary indicates that we could allocate different power to the transmitting 
antennas to achieve more reliable transmission, while still keeping the total energy. Similar 
to the case for Vi(n), one can apply exactly the same algorithm to achieve the ML decoding 
for V2{n). It can be easily seen that these two algorithms have the same complexity. 



2.3 Construction 3 



We will take further efforts to explore the subsets of SU(2). In the following we are going 
to describe a series of unitary constellation V^{n). For the sample program implemented to 
construct V^n), we refer to 0. Now for given integers n > and < k < n, we define 



N = l,N k 



71 



arcsm 



(n— k)i\ 



k = 1,2,- • • ,n, 



M, 



71 



arcsm ■ 



(n — k)n 



2n 



k = 0, 1, • • • ,n - l,M n 



C k (n) = 
Theorem 2.5. Let 



a>k,j bk,i 
—bk,i Ofcj 



\ a k, 



cos • 



(n-k)7iM 



2n 



e N * , h, = sin 



(n-fc)TT i»5 



2n 



-e k 



U CfcW 



fc=0 
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Vz(n) is a fully diverse constellation of^^M k N k elements with diversity product 



k=0 



IT 



nV 3 (n) = sm — 



Proof. Pick two distinct elements A, B e V, 

a b 

A = 



B 



c d 
—d c 



x -b a r 
One can verify that 

| det(A -B)\= det(A - B) = \a - c\ 2 + \b - d\ 2 . 
So if \a\ 7^ |c|, then we have 

\det(A-B)\ = \a-c\ 2 + \b - d\ 2 > (\a\ - \c\) 2 + (|6| - \d\) 2 > 2-2cos^, 

Zi lb 

and it can be verified that the equality holds if there is a k G {0, 1, • • • ,n} such that 
A E Ck(n) and B e Cfc+i(n) or alternatively S e Cfe(n) and A e Cfc+i(n). In the case that 
\a\ = |c|, we will have 

\det(A-B)\ = \a-c\ 2 + \b - d\ 2 > max{\a - c\ 2 , \b - d\ 2 } > 2-2cos— . 



2n 



Therefore for all the cases, we have 



7T 



\det(A-B)\ > 2-2cos— . 

Zilh 

One checks that the lower bound for each case can be reached. So for constellation V 3 (n), it 
follows that 

n vs ( n ) = o( 2 - 2c ° s ^ = sin7r 



2n' 



An 



□ 



Corollary 2.6. When n — > oo, Vz{n) will have have 0(n 3 ) elements and have the diversity 
product O(-). 



Proof. For k — 1, 2, • • • , n, 

N k = 



IT 



arcsm • 



(n — k)i\ 



< 



7T 



Similarly for k — 0, 1, • • • , n — 1, 



M, 



7T 



arcsm • 



(n — fc)7i 



< 



[arcsin sin 



7T 



[arcsin sin 



< 4n. 



< 4n. 
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Hence we derive an asymptotic upper bound for the cardinality of Vz{n), 



|V 3 (n)| = \ C ^ n )\ = J2 MkNk ^ 2n + 16 ( n ~ ^ ^ 

In the following we are going to derive an asymptotic lower bound for |V 3 (n)|. First pick 
two real numbers a, (3 such that 

< a < (3 < 1. 

For k such that an < k < j3n (such k always exists when n is large enough), we have 

N k = > 



sin -j- 

arcsm — 



(n — k)n 
COS - — r, — — 

2n 



sm -p- 

arcsm — 



(1 — Ol)t7 
COS J g-^— 



Utilizing the fact that 

when x is close to 0. When n is large enough, one derives 

(1 — oi)ti 

N k > 4 cos —n. 

2 

Similarly 

M fc > 4 sin - — —^—n. 
Therefore we have an asymptotic lower bound for the cardinality of V^n\ 

|V 3 (n)| > M k* N k> 16(/9 - a) sin ^ ~ ^ cos ^ ~ "^ n 3 > 0(n 3 ). 

on<fc</3n 

Combining the upper bound and lower bound, we conclude that V^{n) has 0(n 3 ) elements 
with diversity product sin(7r/(4n)) = 0(l/n). □ 

This corollary indicates that asymptotically Vz{n) will have much better diversity product 
compared to 0(n), Vi(n), V2(n), because the other three constellations asymptotically will 
have 0{n 2 ) elements and have diversity product O(-). Observe that SU(2) in fact can be 
viewed as 3 dimensional unit sphere. Finding a constellation with the optimal diversity 
product can be converted to a sphere packing problem on this unit sphere. For n points on a 
3 dimensional unit sphere, the largest minimum distance one can hope for is asymptotically 
0(l/n 1 / 3 ). Therefore when n becomes large, asymptotically V$(n) is the best constellation 
over all the subsets of SU(2). Tabled shows the diversity product of constellation Vz{n). 

The ML decoding of V^{n) is also decomposable and it roughly requires 4iV + 2n complex 
multiplications and AN + 2n complex additions, therefore asymptotically it requires O(N) + 
(^(L 1 / 3 ) complex multiplications and additions. One can see from Table^that for reasonably 
small values of n the constellation size grows rapidly. For example when the constellation 
size L is already huge (10770), n is still reasonably small (13). So even if we are dealing with 
a high rate constellation, the decoding process of Vz{n) is still very simple. 
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n 


size 


diversity product 


3 


124 


0.259 


4 


293 


0.195 


5 


582 


0.156 


6 


974 


0.131 


7 


1640 


0.112 


8 


2438 


0.098 


9 


3510 


0.087 


10 


4898 


0.078 


11 


6516 


0.071 


12 


8433 


0.065 


13 


10770 


0.060 



Table 1: the diversity product of V$(n) 

We compare different codes from GPSK signals with transmission rate around 4.5 in 
Figure n One can see that Vs(5) outperforms other constellations even with the highest 
transmission rate. Of course the decoding of Vs(5) is a little more complex than other 
constellations, however the sacrifice in decoding efficiency is worthwhile for the remarkably 
gained performance. 

We compare the GPSK constellations with the Cayley codes jl] in Figure El One can see 
that Vi(44) with transmission rate R = 5.9594 has already a gain of about 2db compared to 
the Cayley code, and the performances of Vs(9) with R = 5.8886 and V^IO) with R = 6.1290 
are even more remarkable. Note that all the GPSK constellations admit the presented 
simple decoding alrogithms. Vi(44) or V2(64) only needs about 8 complex multiplications 
and 8 complex additions to decode one codeword, Va(9) or V^IO) needs about 28 complex 
multiplications and 28 complex additions to decode. 

3 GPSK Constellations from the Real Orthogonal De- 
signs 

In jH] space-time block coding has been discussed for coherent channels. The authors classify 
all the real orthogonal designs and show that orthogonal design codes only exist for 2,4 or 
8 dimension. We show that the normalized version of the real orthogonal codes can be used 
in unitary space time modulation for non-coherent scenarios too. We are going to present 
explicit implementations of the real orthogonal designs, which feature simple decoding al- 
gorithms. For simplicity we describe our idea using a 4 dimensional real orthogonal design: 

/ «i s 2 s 3 

„ _ S 2 Si -s 4 
— S3 S4 Si 

\ — s 4 — s 3 s 2 



S 4 \ 

S3 
S 2 
S! I 



(3.r 
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f V 3 (5)=4.5924 
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10 12 14 16 18 



Figure 1: Performance of different codes from GPSK signals, M=2, N=12 



One natural idea is to transform a 2 dimensional complex orthogonal code 0(n) to a 4 
dimensional real orthogonal code. Recall that 0(n) is a set of matrices 



0{a,b) 



1 

V2 



a 

-b* 



with a, b ranging over ■d = {1, e 27ri//n , • • • , e 27r H n - 1 )/ n }. Let si = Re(a),s 2 = Jm(a),s 3 = 
Re(b), s 4 = Im(b), we will have a 4 dimensional real orthogonal code with the same diversity 
product as 0(n). 

This code also features a simple decoding algorithm. The simple decoding of real orthog- 
onal design codes has been explained in [H], in the following we will discuss the decoding 
process in more details such that it can be generalized to our constellations straightforwardly. 
First we introduce some notations to be used. The superscript T will be the notation for 
a matrix transpose. For two vectors F and G, let F ■ G denote the dot product of F and 
G. For a vector F, Fj denotes the j-th element of F. For a matrix F, Fj denotes the j-th 
column vector. Let B be the mapping 



B(x 1 ,x 2 ,x 3 ,x 4 ) 



( X\ 

X'2 
^3 

\ x 4 



X'2 
—X\ 
— £4 

x 3 



£3 £4 \ 

£4 — £3 

— £1 £2 

-x 2 —xi ) 
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Figure 2: Performances of GPSK constellations and a Cayley code, M=2, N=2 



immediately one checks that 



S 



\ X 4 / 



B(x 1 ,x 2 ,x s ,x 4 ) 



( si \ 

S3 

\s 4 J 



(3.2) 



Assume that the differential unitary space time modulation is used for a wireless commu- 
nication system with M transmitting antennas and N receiving antennas. Let X, Y denote 
the received M x N matrices at time block r — 1 and r, respectively. The ML decoder will 
make the following estimation: 

S = argmin \\Y — SX\\ 2 . 
s " 

We separate X, Y into the real parts and the imaginary parts with 

X = F + iG, Y = P + iQ, 

then we have 



S = argmin (||P — SF\\ 2 + \\Q — SG\\ 2 ) = argrmn 



N N 

0=1 j'=i 
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Utilize Equation ()3.2|) . we have 



S 



arg mm 

S1,S2,S3,S4 



N 

E 

i=i 



N 

+ E 



Qj - B(Gj) 



( si \ 

S3 
\S 4 J 



I 



Since B(Fj),B(Gj) are orthogonal matrices, simple algebraic manipulations will lead the 
above evaluation to 



/ 



arg mm 

Sl,S2,S3,S4 



N 1 

El4 



i=i 



B T (Fj)Pj - \F 



( S! \ 

S2 
S3 
\ S 4 ) 



N 

+E 



Gil 



v 



arg max (si, s 2 , s 3 , s 4 ) • (V* B T (F j )P j 

Si,S2,S3,S4 * ' 

i=i 



v 

E 



|Gi| 



iV 



S3 
\S A J 



B 1 



arg max 



i=i 



where C/j = EE =1 B T (Fj)Pj + EE =1 i? T (Gj)Qj. Since the construction of our code is based 
on 0{n) as described above, one can check that (si,s 2 ) in fact is independent of (s 3 , s 4 ) 
Rewrite si + is 2 



_^ e 2kni/n and S;j + 



5 = arg max ((s x , s 2 ) • (U u U 2 ) T + (s 3 , s 4 ) • (J7 3 , £/ 4 ) T ) 



Sl,S2,S3,S4 



we conclude that the ML decoding of this GPSK constellation is decomposable and can be 
boiled down to the following simple form: 



k 



2vr arg(E/i + iU 2 ) 



n 



I 



2vr arg(f/ 3 + iU 4 ) 



n 



Generally speaking the use of the proposed codes above for a wireless communication system 
with M transmitting antennas and iV receiving antennas will take 8M 2 N real multiplications 
and 8M 2 N real additions to decode one codeword, which is very simple. 

Apparently we can apply the same idea for Vi(n) (i = 1,2,3) to construct GPSK real 
orthogonal constellations. The ML decoder of the corresponding codes will take similar ap- 
proaches as in the 0(n) case, except that one should notice that (si, s 2 ) is not independent 
oof (s 3 , s 4 ) anymore. In this case one can restrict the searching area to be the subsets Cj(n) 
"locally" and apply the similar techniques as we did in the complex orthogonal design case 
to achieve the ML decoding, then the corresponding codes will admit simple decoding algo- 
rithms too. Also we conclude that codes from Vi(n) or V 2 (n) are of the same complexity as 
codes from 0{n). For the codes from V 3 (n), the complexity is increased for high transmission 
rates, however they will have more pronounced performances. 

In the sequel we are going to present a series of 8 dimensional GPSK real orthogonal 
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constellation V^n). Consider the 8 dimensional orthogonal design: 
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se 


s 5 


-S4 
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Sl 
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(3.3) 



Similar to the 4 dimensional real orthogonal constellation case, one can obtain 8 dimensional 
unitary codes by transforming V(n) x W(n), where V or W denotes any one of O, Vi, V 2 , V3 
and x denotes the Cartesian product. For this implementation, Sj can be assigned to be 
the scaled version of the real or imaginary part of "a" or "6" as in 4 dimensional case. 
However as in the complex orthogonal case Vs(n) motivates us to explore more densely 
packed constellations. 

The problem of constructing an 8 dimensional real orthogonal code with the maximal 
diversity product is equivalent to the packing problem on a 7 dimensional unit sphere. There- 
fore any currently existing results in packing problem on a 7 dimensional unit sphere can be 
"borrowed" for the real orthogonal code construction. However decoding such codes using 
exhaustive search will be unpractical for high transmission rates. In the sequel we are going 
to present a series of 8 dimensional GPSK orthogonal codes V 4 (n) featuring simple decoding 
algorithms. 

The basic idea is that instead of considering the problem of packing the real vector 
s = (si, s 2 , • • • , s 8 ) on a unit sphere, we consider the problem of packing the complex vector 
z = (zi, z 2 , z 3 , Z4) with constraints Y^j=i \ z j\ = 1> where Zj = s 2 j-i + is 2 j. Since the distance 
of two complex vectors satisfies 

4 4 

3=1 j=i 

we can pack the amplitude vector (\z±\, \z 2 \, \z 3 \, |z 4 |) first, then pack the argument vector 
(arg(zi), arg(z 2 ), arg(z3), arg(z 4 )). To explain this implementation in more details, we repre- 
sent z using the polar coordinates: 

zi = cos#ie ni , 

z 2 = sin9 1 cos9 2 e n2 , 

z 3 = sin 61 sin 9 2 cos 9 3 e n3 , 

z^ = sin 6 1 sin 8 2 sin # 3 e i74 . 

Let 9j run over rrij + 1 evenly distributed discrete values from to 71/ 2: 

and let jj run over rij evenly distributed discrete values from to 2n: 

7j G {k — I A; = 0, 1, • • • , rij — 1}. 

rij 
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Now we have a finite set of complex vectors. 

Let v = (v%, v 2, v 3 , V4) and w = (iy 1; w 2 , w 3 , w 4 ) be two distinct resulting complex vectors. 
If \vx I 7^ |wi|, one can check that 

I v — w 1 2 > 4 sin 2 (it / {Ami ) ) ■ 

If \v\\ = \w\\ but |f2| 7^ \w 2 \, we wm have 

|v - w| 2 > 4(1 - |^i | 2 ) sin 2 (7r/(4m 2 )). 

In the case that \v\ \ = \w\ \ and |t»2 1 = 1^2 1 but | ^3 1 7^ \ws\, similar algebraic calculations will 
give a lower bound 

|v — w| 2 > 4(1 — |fi| 2 — |f2| 2 ) sin 2 (7r/ (Arris)). 

We further consider the case \vj\ = \wj\ for j = 1, 2, 3, 4. In this case, if Vj 7^ Wj for some j, 
we can have 

|v — w| 2 > A\vj\ 2 sin 2 (7r /rij). 

It is easy to check that the lower bound can be reached for some special pair of z, w in all the 
cases. Based on the 8 dimensional real orthogonal design (jH2J), assign s^j-i = Re(zj) and 
S2j = Im(zj) for j = 1,2,3,4. Now we have a finite unitary constellation whose diversity 

product is 

min a/|v — w|/2. 

v,w 

However if we choose mj,rij arbitrarily, often we will end up having a constellation with 
a small diversity. In the sequel we are going to add suitable constraints to the choices of rrij 
and rij to guarantee that we have a large diversity product. We present the algorithm step 
by step. For the sample program of this algorithm, we refer to 

Algorithm to Construct V^n): (given the input n) 

1. Fix mi = n and let Q\ run over n + 1 even distributed discrete values from to 71/ 2 
(consequently \z±\ runs over n + 1 discrete values from to 1). 

2. For any fixed \zi\, take m 2 to be the largest integer such that 

4(1 - |^i| 2 ) sin 2 (7r/(4m 2 )) > 4 sin 2 (7r/(4n)). 
Let 62 run over m 2 + 1 evenly distributed discrete values from to tt/2. 

3. For any fixed \zi\, \z2\, take m 3 to be the largest integer such that 

4(1 - |^i| 2 - |z 2 | 2 )sin 2 (7r/(4m 3 )) > 4 sin 2 (7r/(4n)). 
Let 83 run over m 3 + 1 evenly distributed discrete values from to tt/2. 

4. For any fixed \zi\, \z2\, \zs\, take rij to be the largest integer such that 

41^- 1 2 sin 2 (tt/t^) > 4sin 2 (7r/(4n)). 
Let jj run over rij evenly distributed discrete values from to 2tt. 
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Comparison of Different Dimensional Constellation 

10° F 1 1 1 1 
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SNR(db) 

Figure 3: Performances of a GPSK constellation and a code with optimal diversity product, 
N=2 

5. The steps above result in a finite set of complex vectors z = (z%, z%, z%, 24). Based on the 
8 dimensional real orthogonal design (|3.3jh assign S27-1 = R e ( z j) an d S2j = Imizj) 
with j = 1,2,3,4. Now we have a fully diverse unitary constellation with diversity 
product sin(7r/(4n)). 

We call this constellation Vjji). Recall that for the packing problem with n points on 
a 7 dimensional unit sphere, asymptotically the largest minimum distance one can hope 
for is 0(l/n 7 ). The following theorem indicates that asymptotically V&{n) are the best 
constellations from 8 dimensional real orthogonal designs. Although the proof of this theorem 
is quite tedious, the basic idea is very similar to Corollary 12.61 thus we skip the proof. 

Theorem 3.1. Va(ti) is a fully diverse constellation with diversity product sin(7r/(4n)). 
When n — > 00, V&{n) will have have 0(n 7 ) elements and have the diversity product Oil/n). 

Applying the same analysis for the constellation V${n), we conclude V±{n) is also decom- 
posable and the decoding complexity for one codeword will be O(N) + 0(L 3 ^ 7 ). In FigureEl 
we compare the performance of 14(2) of 103 elements with a 2 dimensional constellation of 
3 elements: 

{h,AB}, 

where A = diag (e l27r//3 , e l27r//3 ) and B = A 2 . It is well known that this 2 dimensional con- 
stellation has the optimal diversity product over all the 2 dimensional unitary constellations 
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with 3 elements. Since a large number of transmitting antennas guarantee that the full diver- 
sity at the transmitter side can be utilized more efficiently (see Inequality II. 2|) . it is not too 
surprising to see that the 8 dimensional GPSK constellation performs better at high SNR 
region. The figure shows that from 6db, 8 dimensional constellation begins to outperform 
the 2 dimensional one. Around Block Error Rate (BLER) of 10~ 6 , the performance gain is 
about 5db, which is very pronounced considering the decoding of V±(2) is relatively easy. 

4 Conclusions and Future Work 

The complex and real orthogonal coding schemes admit simple decoding algorithms. Based 
on these schemes, we generalize one dimensional PSK signals and explicitly construct GPSK 
unitary space time constellations. These constellations can be viewed as higher dimensional 
generalizations of one dimensional PSK signals and theoretical analysis shows that their 
decoding procedures are decomposable, i.e., the demodulation of these codes can be boiled 
down to one dimensional PSK demodulation. Therefore our constellations have very simple 
decoding procedures. For some of the resulting codes (for example, Vi(n), V2(n)), the com- 
plexity of ML decoding does not even depend on the transmission rate. We reallocate the 
power among the antennas (meanwhile keeping the total energy) to optimize the diversity 
product. Numerical experiments show that our codes perform better than some of the cur- 
rently existing comparable ones. For the sample programs regarding how to construct the 
proposed constellations, we refer to [3]. 

The theoretical analysis JS] shows that the complex and real orthogonal design constel- 
lations only exist on 2, 4 or 8 dimension. Based on these schemes we only propose 2, 4 and 
8 dimensional GPSK constellations. We are trying to generalize the idea of decomposable 
GPSK constellations for any dimension. Constellations with simple decoding algorithms 
based on the generalized orthogonal designs [7j are under investigations. 
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